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We investigate entanglement distribution in pure-state quantum networks. We consider the case 
when non-maximally entangled two-qubit pure states are shared by neighboring nodes of the net- 
work. For a given pair of nodes, we investigate how to generate the maximal entanglement between 
them by performing local measurements, assisted by classical communication, on the other nodes. 
We find optimal measurement protocols for both small and large ID networks. Quite surprisingly, 
we prove that Bell measurements are not always the optimal ones to perform in such networks. 
We generalize then the results to simple small 2D networks, finding again counter-intuitive optimal 
measurement strategies. Finally, we consider large networks with hierarchical lattice geometries 
and 2D networks. We prove that perfect entanglement can be established on large distances with 
probability one in a finite number of steps, provided the initial entanglement shared by neighboring 
nodes is large enough. We discuss also various protocols of entanglement distribution in 2D networks 
employing classical and quantum percolation strategies. 

PACS numbers: 03.67.-a, 03.67.Bg 



I. INTRODUCTION 

Quantum Networks play a key role in quantum in- 
formation processing. In such networks, quantum states 
can be prepared initially and shared between neighboring 
nodes (or stations), i.e. entanglement can be generated, 
and this resource is then to be used for quantum com- 
munication [3, 3], or distributed quantum computation 
Q involving arbitrary nodes of the network. One of the 
main tasks is then to design protocols that use the avail- 
able quantum correlations to entangle two nodes of the 
network, and to optimize these protocols in terms of final 
entanglement and probability of success. 

A set of quantum repeater stations, for instance (see 
Fig. [T^), can be considered as a ID quantum network, 
where the aim is to establish quantum communication 
over large distances [QT9i| . It is well known that the sim- 
ple entanglement swapping procedure can achieve 
this goal, but (except for the unrealistic case of per- 
fect resources and operations) the probability of ob- 
taining entanglement between the end-nodes of such a 
network decays exponentially with the number of re- 
peaters. This problem can be overcome by the more 
sophisticated quantum repeaters protocols which 
intersperse "connection steps" (entanglement swapping) 
with purification steps and require only polynomial de- 
cay, thus opening the way for feasible long-distance quan- 
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tum communication. 

In two-dimensional or higher dimensional lattices of 
large size, a perfect connection between any two nodes is 
possible with a probability that is strictly greater than 
zero, even with imperfect resources. This can be achieved 
by the so-called classical and quantum percolation strate- 
gies [ll|, in which initially, or after some preparatory 
measurements, respectively, one converts all bonds into 
singlets with a probability p. This result is very encour- 
aging, but remains of little use for finite and small net- 
works [12, The aim of the present paper is twofold: 
first, we investigate and derive optimal local measure- 
ment protocols for simple networks of finite size. In par- 
ticular, we consider certain ID and 2D networks of nodes 
that consist of z qubits, where z is the number of neigh- 
bours. Neighboring nodes share partially entangled pure 
states. We apply then local quantum operations to the 
nodes, assuming that these operations are noiseless. We 
first address the question of optimal entanglement prop- 
agation, or in another words entanglement swapping, in 
small networks consisting of three or four nodes only. The 
insights obtained for these simple situations are then used 
as building blocks for larger ID and 2D quantum net- 
works, as well as networks with hierarchical geometry. 

Our second aim is to discuss examples of hierarchical 
"diamond" and "tree" lattices in which perfect entan- 
glement on arbitrary large distances can be achieved 
in a finite number of steps (measurements). Provided 
that sufficiently large but not necessarily maximal 
entanglement is present, we can convert connections 
along a given line into prefect singlets. Finally, we 
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consider various kinds of percolation strategies: the one 
presented in Ref. [ij], which employs a change in the 
lattice connectivity due to quantum measurements, and 
a different one in a triangular lattice, where the optimal 
singlet conversion strategy is used. Both of these 
protocols essentially demonstrate that the quantum 
percolation thresholds are lower than their classical 
counterparts. Equally interesting we propose to use 
the optimal singlet conversion strategy to transform a 
square lattice into two independent square lattices of 
doubled size, for which the percolation probability is 
larger than in the original lattice. 
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FIG. 1: Notation and examples of ID networks: (a) the stan- 
dard quantum repeater scenario; (b) entanglement swapping; 
(c) a two-repeater system. 



Outline In Section II we fix the notation and define 
the figures of merit used for evaluating the measurements 
efficiency: the concurrence (C), the so-called worst case 
entanglement (WCE) and the singlet conversion proba- 
bility (SCP). In Sec. IIIII we describe the strategies max- 
imizing these quantities for some ID networks, starting 
from a simple one-repeater configuration, consisting of 
two bonds with two imperfectly entangled pairs on them. 
Interestingly, there exists a strategy that conserves the 
averaged singlet conversion probability [l^ ; the protocols 
however does not scale with the number of repeaters, as 
expected. The second subsection of Sec. lIIII deals with the 
problem of two repeaters, that is three bonds. Here the 
optimization of the SCP is much more complex: in some 
conditions we obtain that the optimal measurements do 
not correspond to a Bell measurement, i.e. von Neumann 
measurements in a Bell basis of orthonormal maximally 
entangled states. This result is somewhat an alog ue to 
the recent result by Modlawska and Grudka [1^|, who 
have demonstrated that non-maximally entangled states 
can be better for the realization of multiple linear optical 
teleportation in the scheme of Knill, Laflamme and Mil- 
burn The last part of this section deals with large 
ID network {i.e. in the limit of infinite size network). 
Here we prove that the probability of establishing entan- 
glement over large distances decays exponentially. We 
present optimal strategies for the concurrence and the 
WCE, and upper bounds for the SCP. 

In Sec. IIVI we turn to the simplest small network in 
2D: a square. We obtain similar results as in the case of 
two repeaters in ID, indicating that Bell measurements 
not always provide the best protocol. In Sec.|V]we apply 
the results of previous sections to networks of large size 
and hierarchical geometry, that is, lattices that iterate 
certain geometric structures, so that at each level of iter- 
ations the number of nodes, or the number of neighbors 
changes. We consider two kinds of hierarchical lattices: 
first we discuss the so-called "diamond" lattice, for which 
we prove that for sufficiently large initial entanglement, 
one can establish perfect entanglement on large scales 
{i.e. some lower levels of iteration) in finite number of 
steps. A somewhat simpler result holds for the simplest 
possible double Cayley tree lattice, in which in each step 
of iteration each bond branches into two. For such lat- 
tices, if the initial entanglement is large enough, perfect 



entanglement can be established at each level of iteration. 

Finally, in Sec. I VII we consider genuine 2D lattices. 
First, using similar method as in Sec.|V]we show that for 
a sufficiently broad strip of a square lattice, we can con- 
vert connections of a given line along the strip into a line 
of perfect singlets, provided, of course, that initially an 
imperfect, but sufficiently large entanglement is present. 
Second we reconsider percolation strategies and discuss 
the example of hexagonal lattice with double bonds from 
Ref. and a triangular lattice with variable bonds. In 
the first of these examples quantum measurements lead 
to local reduction of the SCP, but change the geometry 
of the lattice, increasing its connectivity and thus the 
classical percolation threshold. In the second example 
we use a protocol optimizing the SCP to transform the 
original lattice to a new one with the same geometry, 
but with a higher probability p of getting a singlet on a 
bond. Similarly, we discuss a different type of strategy, 
where by using the optimal singlet conversion protocol we 
transform a square lattice into two independent square 
lattices with the same mean SCP as the initial one. We 
prove that the classical probability of connecting a pair 
of neighboring points in the initial lattice (two neighbor- 
ing point from the two lattices) to another such pair is 
strictly larger for the case of two lattices. We conclude 
then in Sec. |VlIl 



II. PRELIMINARIES: NOTATION AND BASIC 
NOTIONS 

A pure state of two qubits is represented by a solid line 
in the figures and is written (except when specified) as 

|a) = V^|00) + V^|11), (1) 

where ao + ai — 1 and ao > ai (it is assumed that 
local basis rotations are performed whenever necessary 
for the states to be written in that way). This defines 
the Schmidt decomposition of a pure state of two qubits, 
while ao and ai are their Schmidt coefficients. 
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1. Entanglement Swapping 

A basic operation for propagating entanglement over 
larger distances is the so called "entanglement swap- 
ping" , see Fig. [T]d. A POVM is described by n positive 
operators G M4(C) satisfying the completeness rela- 
tion Yl=i E,n = I4. Let 

Pm = tiBC ((12 <^ E„^ (E) I2) \af3){af3\ ) 

be the resulting state of the measurement M , which oc- 
curs with a probability Pm = tr(pm). We consider in 
that paper projective measurements only, i.e. n = A and 
Em ~ \'Um){um\ for somc normalized state For 
such measurements the smallest Schmidt coefficient of 
Pm is A™ = min{eig(p„)}/p„, where Pm = tr^(p„i), or 
equivalently 



(2) 



Considering the following map from to A/(C, 2) 
1 



a 



hj=0 



aio an,' 



one can show that pm is now equal to X^X^, with Xm — 
aum /?■ 

The concurrence of a state f is by definition C'{ip) = 
2|det(^)|. Therefore, the concurrence Cm, the smallest 
Schmidt coefficient Xm and the outcome probability pm 
are explicitly given by 



C„ 



— '~^\Um), (,4aj 



1 

Pm — ^ ^ ^iPj |^m,zj| ■ 



(4b) 
(4c) 



We now turn to Bell measurements. Starting from the 
computational basis { |0), |1)} of a single qubit, we define 
the new ortho normal basis { |t), li)} 



and the Bell vectors 

|5,±)= ITT>± IU> 



V2 



and l^t^^ 



U eU{2), 



ITi)± lit) 



V2 



(5) 
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Two specific bases play a key role in this paper: the com- 
putational or "ZZ" basis, where the vectors |t) and \i) 
for both qubits are the eigenvectors of the Pauli matrix 
CTz, and the "XZ" basis, where the first basis is chosen 
as being the eigenvectors of ax. Although we could in 



principle parameterize the Bell states in that way, calcu- 
lations are much easier and clearer in the "magic basis" 
defined as [IT] 



($1, $2, $3, $4) = (I2 



(7) 



which is nothing but the usual Bell basis with some spe- 
cific phases. In this basis, the concurrence of a state 
Im> = ELiAiil^*) simply reads C{p) = \J2t=it^l\- I* 
follows that the coefficients fii of a Bell state (whose con- 
currence is 1 by definition) have all the same phase; hence 
we can choose them as being real. Let a set of four such 
states {pm}, so that the matrix {pm,i) belongs to 50(4). 
Then the probabilities given in Eq. (|4c| read 



Pm Pmin (^Pm,l Pm.'l) ^" Pmax (^Pm,3 + PmA) , (8) 

with 

Pmin = ^ and Pmax = ^ • (9) 

We emphasize the fact that (given two states a and (3) , 
the outcome probabilities completely characterize a Bell 
measurement, since Xm depends only onp^ for C{um) = 
1, see Eq. g]). 

2. Figures of Merit 

We describe here three figures of merit used to evalu- 
ate the usefulness of an entanglement swapping protocol: 
the concurrence, the singlet conversion probability 
(SCP) and the worst-case entanglement (WCE). All 
these figures of merit take value in the interval [0, 1]. 

Concurrence The average concurrence of a measure- 
ment M is defined as 

CM=Y,PmCm, (10) 

m 

where Cm is the concurrence of the outcome m. 

WCE The idea of the WCE is to find a measurement 
optimizing the entanglement for all its outcomes. Tak- 
ing the smallest Schmidt coefficient as the entanglement 
measure we define the WCE as 



Wm = 2min{A„i}. 



(11) 



SCP We consider here the probability of conversion 
of a given state into a perfect singlet. A result of ma- 
jorization theory [l^, [l^ tells us that a state |a) ([T|) 
can be converted into a singlet by LOCC with maximal 
probability 2q;i; this is the "Procrustean method" of en- 
tanglement concentration described in [13] ■ We define 
the average SCP for a measurement M as 



S 



M 



A. 



(12) 



4 



where Am is the smallest Schmidt coefhcient of the out- 
come TO. Since this figure of merit is used for different 
systems, we sometimes use the following notation for clar- 
ity: 

cWf \ 

where N means the number of repeaters of a ID chain 
consisting of iV -|- 1 states ai, a2, ■ • ■ , ctN+i, as depicted 
in Fig. [T^. 



The result of the maximization of the WCE is summa- 
rized in the following result: 

Result 2 Best WCE strategy for one repeater 

The maximum value of W for a one-repeater 
system is reached by the Bell measurement in 
the XZ basis, with 

W„.ax = Wxz = 1 - Vl- 16aoai/3o/3i- (14) 



III. ID NETWORKS 

Before studying the two-dimensional networks, it is 
worth looking at systems made of one or two repeaters 
only. In fact, some interesting properties of these small 
systems can then be used in more elaborated strategies 
for larger networks. For instance, the important fact 
that the SCP does not decrease after one measurement 
f jjlll A3P allows one to get better results for the perco- 
lation on honeycomb lattices [T]| . Another important 
and surprising result is that Bell measurements are not, 
in general, the measurements that maximize the SCP 
( mil B 2[) . although they are the best ones for the aver- 
age concurrence and the WCE. Previous results on ID 
networks can also be found in Rcfs. (2ll. [23. 



A. One Repeater 

We consider in this section a system consisting of two 
states a and /3 joined by a single repeater, see Fig. [TJd. 
We first prove a general statement on Bell measurements, 
and then describe the measurements that maximize our 
three figures of merit. 



1. Bell Measurements and Outcome Probabilities 

The following result is very useful when trying to 
maximize the SCP over the set of Bell measurements 
(the proof is given in App. [X| : 

Result 1 Outcome probabilities for a one-repeater Bell 
measurement 

Let {xrn} be four real numbers that add up to 
one and that lie in the interval [pmin, Pmax]- 
Then there exists a Bell measurement whose 
outcome probabilities pm are equal to Xm ■ 



2. Maximizing the Concurrence and the WCE 

It is clear from Eqs. (|4al [T0|) that any Bell measure- 
ment, i.e. C{um) = 1 Vto, maximizes the average con- 
currence, and therefore 



2^/ aoai/3o/3i. 



Proof. (By contradiction). The Bell states \um) 
XZ basis are given by the columns of the matrix 



the 



XZ 




hence Pm = 1/4 and 2Am = 1 — \/T^^T6aoai^3o^ Vto. 
Now suppose that there exists a measurement M de- 
scribed by the set {E„i = |um)(um| with n > 4, 
such that Wm > Wxz- Then each Am has to be strictly 
greater than the smallest Schmidt coefficient of the out- 
comes in the XZ basis. Thus, from Eq. ^ 



det{pm) > pli-iaoaiPoPi Vto. 



(15) 



Since det(pm) = aoai/So/^i | det(tim)p, the summation 
over TO of the square root of Eq. yields 



J2 |det(um)| > 2. 



(16) 



m— 1 



But the concurrence of a (normalized) state is smaller 



or equal than one, hence 2|det(u„j)| < 



More- 



over, taking the trace of the completeness relation for 
the operators impHes J2m=i = 4. There- 

fore J2m=i I det(um)| < 2, which is in contradiction with 
Eq. (fTB]) and concludes the proof. □ 



3. Maximizing the SCP 

The following result gives the maximum value of the 
SCP for one entanglement swapping step: 

Result 3 Best SCP strategy for one repeater 

The measurement that maximizes S for a 
one-repeater configuration is the Bell mea- 
surement in the ZZ basis, and 



Szz = 2 min{Q;i,/?i}. 



(17) 



(13) 



Proof. Two kinds of outcomes appear when performing a 
Bell measurement in the computational basis: two of the 
outcome probabilities are equal to Pmax, while the other 
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two are equal to Pmin- Putting these values into Eq. (|4|) 
one finds the corresponding smallest Schmidt coefRcients: 

w X w N mm{aof3i, ai(3o} 

A(Pmax) = 7^ , A(pmin) = ^ , (18) 

whence Szz — 2 min{ai,/3i}. Consider now that we are 
allowed to perform some arbitrary unitary not only on 
BC, but on ABC. We are in presence of a bipartite sys- 
tem, and results of majorization theory apply: the SCP 
of this system is at most 2 /3i . A similar construction for 
qubits B, C and D tells us that the SCP is at most 2q!i, 
so that the final SCP cannot exceed twice the minimum 
of ai and /3i. □ 
Remark Setting a = /3, one sees that the SCP does 
not decrease after one entanglement swapping; this is the 
"conserved entanglement" described in [14| . 



B. Two Repeaters 

We consider a system of three states on which we per- 
form two consecutive entanglement swappings, as shown 
in Fig. [!}:, and we describe the measurements that max- 
imize the three figures of merit. 



1, Maximizing the Concurrence and the WCE 

The maximization of these two figures of merit is triv- 
ial for a two-repeater configuration once one knows the 
results for the one-repeater system. First, any Bell mea- 
surement maximizes the average concurrence of the re- 
sults of the two measurements. This will be generalized 
and proved for any number of repeaters in ^III C II Then, 
in order to maximize the WCE, we simply have to per- 
form XZ measurements at both repeaters. In fact, if 
we perform any other measurement on the first repeater, 
then at least one resulting state will be less entangled 
than the XZ results, and this reflects on the WCE of the 
second measurement (which has to be a Bell measure- 
ment in the XZ basis from Result [2]) . 



2. Maximizing the SCP 

After the first measurement, we get four resulting 
states \ipm) with probabilities Pm- From Result [3] we 
know that for any outcome, the second measurement has 
to be done in the ZZ basis. Hence, we have to find the 
first measurement M that maximizes 



Kv) 



'S'm (ao, ^0, 7o) = 2 ^ 



Pn 



^WmS, 7l}- (19) 



We first maximize this quantity over the set of Bell 
measurements (which, as we will see, leads to the best 
strategy for a large range of entangled states a, /3 and 7), 
and then we present some numerical results showing that 
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FIG. 2: Representation of tlie 
min{/(p), g{p)} governing tlie SCP 
ments in a two-repeater configuration. 



function h[p) = 
after Bell measure- 



non-Bell measurements sometimes provide better results. 

Bell Measurements We fix the states a, [3 and 7 and 
we consider the SCP as a function of the outcome prob- 
abilities only: 

S{{Pm}) = ^min{/(p„i), g{p,n)} = ^ /i(Pm), (20) 



where /(p) = 271P and g{p) = p — \/p^ — aQai(3o(3i. One 
can show that g'{p) < and g"{p) > Vp G [pmin, Pmax]- 
A typical plot of h{p) is shown in Fig. [21 and the value 
p* at which the functions / and g cross each other is 



1 apfXiPoPi 



7o7i 



(21) 



It is sufficient to maximize the function over the pos- 
sible probability distributions, since Result [T] insures the 
existence of a Bell measurement leading to this optimal 
distribution; we recall that the probabilities have to be 
chosen in the interval [pmin , Pmax] ■ Let us give two nec- 
essary conditions that have to be satisfied by the best 
probability distribution (they can be proven rigorously, 
but a look at Fig. [5] may be clearer): 

• if the set {pm} maximizes S, then all probabilities 
lie either to the left of p* or to its right. In fact, 
suppose for example that pi -|- 2e < p* < p2 — 2£, 
and choose pi — pi + e and P2 — P2 ~ ^ (with 
< e <C 1 as it should be). The constraints on 
these new probabilities are clearly satisfied if it was 
the case before, and a better SCP has been found. 



• if pi and p2 are such that p* + 
Pmax— 2£, then the choice pi = pi 



2e < pi < p2 < 
-e andp2 = P2+£ 
gives rise to a strictly greater SCP (this comes from 
the convexity of g). 

It is now simple to maximize the SCP of two repeaters, 
and one sees that the value p*, with respect to Pmin 
and Pmax, plays a crucial role in the choice of the best 
probability distribution. In fact, we have to distinguish 
four distinct cases; see results in Tab. HI We notice that 
ZZ measurements lead to the maximum SCP whenever 
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value of p* 


{pm} maximizing S'''^' 


Jr" J-* ITI 1 11 


L/'^mini /-^miii 7 7-^max; /-'max J 


Pmin < P* < (1 — Pmax)/3 


{p*, p*, Pmax, 1 - 2p* — Pmax} 


(1 -pmax)/3 <p* < 1/4 


{P*, P*, P*, 1 - 3p*} 


P* > 1/4 


{1/4, 1/4, 1/4, 1/4} 



TABLE I: Maximization of 5"'^' over Bell measurements, see 
text for details. 



P* < Pmin, while the XZ ones are the best strategy for 
p* > 1/4. So far, we have maximized the SCP for two 
repeaters supposing that the first measurement was to 
be done on the states a and /3. But what happens if we 
start from the right side? It appears that the maximum 
SCP depends, in general, on the order of the mea- 
surements and that performing the first measurement 
where the states are more entangled yields better results. 

General Measurements (Numerical Results) The 
question is to check if some non-Bell measurements yield 
a better SCP than the results of the last paragraph. Since 
the concurrence of the states used for entanglement swap- 
ping can now take any value between and 1, we cannot 
consider 5* as a function of the outcome probabilities only. 
But for a fixed concurrence C < 1 one sees that 



9{C, P) = P- \Jp^ - aQai/3oPi < g{p) Vp. 

Writing the corresponding variables of non-Bell mea- 
surements with a bar, we have that p* < p* and 
g{C, p*) < g{p*)- Therefore, one can check that Bell 
measurements are indeed the best ones, except, possibly, 
when Pmin < P* < (1 — Pmax)/3. The key fact about Bell 
measurements in that case is that we cannot chose three 
outcome probabilities to lie on p* , since the fourth one 
would be greater than Pmax- But the range of possible 
outcome probabilities depends on the concurrence: for 
example, from Eq. (Pc)l and for C(um) = 0, we have that 
Prn G aoPo], or more generally 



Pm G [Pmaxi Pn 



[Pmaxi Pn 



(22) 



Hence, and this is confirmed by numerical results, a bet- 
ter strategy is to perform a measurement such that three 
outcomes probabilities are equal to p* , and that the con- 
currences of the states are the largest ones satisfying 
Pmax = 1 — 3p*. Our numerical evidence shows that 
Bell measurements do not always maximize the SCP, see 
Fig. El 



C. Large ID Chains 

We consider the system of Fig. [T^ that consists of N 
repeaters joining -I- 1 states. For simplicity, we choose 
the states ai as being identical: \ai) = \(p) Vi. We show 



S'mlx(ao,/3o,7o) 
0.4 




FIG. 3: SCP for a system of two repeaters, with /3o = 70 = 
0.7. Numerical results (dashed line) show that there exists a 
better strategy than Bell measurements (solid line) for ao £ 
[ai, 02]. The values ai and a2 are such that p*(ai) = (1 — 
Pmax(ai))/3 and p*(a2) = Pmiii(a2)- 



in this section which strategies yield the optimal solution 
for the concurrence and the WCE, and for the SCP we 
give an upper bound to its maximum value and some 
results for XZ and ZZ measurements. 



1. Maximizing the Concurrence and the WCE 

A direct generalization of Eq. (|^a|) for N repeaters 
yields for the concurrence [1^: 



C(^) = ^ 2 I det (X 



{m.lJ 



(23) 



where X 



(p, and the states \ur, 



are associated with the measurement result rui of the i-th 
repeater. Then the maximization of C*-^' reads 



maxjC^^H = |det(£)|^+i max 



det ($u,„i 



{ E 2""^ 

{r«.} 

)|} = |2det(^)r\ 



(24) 



where $ — l2/'\/2 corresponds to a maximally entangled 
state. For states p which are not perfect singlets, the 
concurrence decreases exponentially with N: 



(25) 



The same arguments as for the systems of one or two 
repeaters hold for the WCE, so that XZ Bell measure- 
ments have to be performed on each repeater in order to 
maximize it. 



2. Maximizing the SCP 
A similar formula as Eq. (|23p for the average SCP is 
^2min{eig(X,„^Xt^J}. (26) 

{717;} 
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CS 



xz 



zz 



N/2 



TABLE II: Asymptotic behavior of the SCP for a ID chain of 
A'' S> 1 repeaters. Three specific measurement protocols are 
studied: conversion of aU states into singlets (CS), XZ and 
ZZ measurements. 




FIG. 4: Operations on a square to obtain an entangled pair 
on the diagonal: first two measurements, then distillation of 
the resulting states a and 



Contrary to the maximization of the concurrence, we can- 
not find here such an easy way to calculate the maximum 
value of S, but we already can say a few words about the 
SCP for a ID chain with a large number of repeaters: 

• Since S is always smaller than or equal to C, it is 
upper bounded by 



(27) 



• After ^ 1 measurements, the entanglement of 
the resulting states is expected to be, in average, 
very small, so that the SCP and the concurrence 
could be related by S" w C^. Hence we may have 
the asymptotic behavior S^^'^ ~ (AipQipi)^ . 

Even if we do not have the protocol that maximizes 
the SCP, we present here three specific strategies, as the 
results are instructive. The first and simplest one con- 
sists in trying to convert each state into a singlet, and 
then to establish a perfect connection between the end- 
qubits of the chain. In the second strategy we perform 
XZ measurements at all stations, and from i jlll A2I we 
know that all resulting states have the same amount of 
entanglement. We indeed find the exponentional decay of 
the SCP related to the one of the concurrence. Finally, 
in App. [B] we derive the explicit formula for ZZ mea- 
surements on a chain of any number of repeaters, which 
yields an decay of the SCP which is quite close to the up- 
per bound given in Eq. (|27p . The asymptotic behaviors 
are summarized in Tab. [ITl 



IV. THE SIMPLEST 2D NETWORK: A SQUARE 

We study in this section a square made of four iden- 
tically entangled states, see Fig. S) This is clearly one 
of the simplest possible 2D networks. The operations we 
perform consist of three steps: a first measurement Mi 



yielding some outcome a, then a measurement M2 de- 
pending on a and giving another state /3, and finally a 
distillation of these two states to get a final state ifj. The 
goal is of course to get ip as entangled as possible, given 
the states 



1. Distillation 

Majorization theory [l^ tells us how entangled the 
state i/j can be. Without loss of generality we choose 
Q!o ^ Po and the majorization criterion reads 



{aoPo, aoPi, ai/3o, ai/3i) -< (V'o, V'l, 0, 0) 



(28) 



whose only non-trivial inequality is ao/^o ^ "00- Since we 
are looking for a state ip that is as entangled as possible, 
we know its greatest Schmidt coefficient: 



ao/3o|- 



(29) 



2. Maximizing the Figures of Merit 

Arguments used for ID networks still hold here, so 
that one has to perform Bell measurements and XZ mea- 
surements to maximize the concurrence and the WCE 
respectively. It is worth pointing out that a perfect sin- 
glet ip can be established with probability one after two 
XZ measurements and a distillation if ip satisfies 



0.65. 



(30) 



Thus we consider that if is less entangled that tp* 
since we already know how to get a singlet for tpo < ipQ. 
We proceed in two steps for maximizing the SCP: we 
first look at the subproblem of maximization over the 
measurements M2 for a given outcome a, and then we 
provide some numerical results for the whole square. 

Second Measurement We first notice that a singlet 
can be obtained by a XZ measurement with probability 

one if ao < = ^1 + ^1 — (4 ipoipi)"^^ . Then, label- 
ing by m the resulting states /? of the measurement M2, 
we can write the function to be maximized as 



'a/ 



2^p, 



1^, ao/3m,o| 

. r ao-ail 
mmUim,!, — > 

L Ian J 

m 

S^Licio) + ao Sf} ((^0, Vo. 77—) : (31) 
V z an J 



= 2a 



so that all results of Sect. IIII Bl can be applied. The 
three quantities Pmim Pmax and p* used in that sec- 
tion are now pmi„ = ipo'fi, Pmax = {ipl + (p\)/2 and 



8 



Hp) 
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5° 
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P 0.25 Pmax 

FIG. 5: Typical plot of the function h{p) governing the SCP 
of the square. 



FIG. 6: SCP for a square made of four states (p. Numerical 
results (dashed line) show that Bell measurements (solid line) 
do not lead to the optimal solution for ipo > if* 0.664. 



value of p* 


{pm} maximizing 


<po ^ I 


{Pmin; Pmini PmciX^ Pmax} 


fo > fo 


{p*, p*, Pmax, 1 - 2p* - Pmax} 


fo < fo < 


{p*, P*, P*, 1 - 3p*} 




{1/4, 1/4, 1/4, 1/4} 



TABLE III: Maximization of S over Bell measurements, see 
text for details. 



p* — ipoLpi Q!o/(q;o — cti)- Since p* is greater than Pmin for 
all states a and ip, it follows that S'^g.x is reached by Bell 
measurements except when p* G ]pmin, (1 — Pmax)/3[. 

First Measurement The function to maximize over 
the measurements Mi is 

^Ml = X^P™ '5'^ax("m,0, fO, <Po)- (32) 
m 

For Bell measurements, since the Schmidt coeffi- 
cient am,o depends on pm only, we can write S^j_^ = 
h{pm)- Here we make a slightly abuse of notation, 
since we again use h{p), as in ijlllB 21 Actually, the shape 
and properties of the function h(p) discussed here and in 
i jlll B 2l are very similar. Therefore, all arguments used in 
that section for the maximization of the SCP apply here, 
too. The plot of h{p) is shown in Fig. [5] The quantity 
that corresponds to p* is now written p* and its value is 

where a\ = 1 — With these definitions, one can check 
that for all greater than (/Jq, we have Pmin ^ P*{'^o) ^ 
1/4 and that p* — > Pmin when i^sq ~> 1, whence the best 
measurements for nearly unentangled states Lp are the 
ZZ ones. As for the system of two repeaters, performing 
Bell measurements is not the best choice when it is not 
possible to get three of the four outcome probabilities to 
be equal to p* (but this possible when 'Pa < ~ 0.664, 
see Fig. ^ . Finally, we summarize the results in Tab, [ml 
and the similarity with Tab. [T] is immediate. 



V. HIERARCHICAL LATTICES 

In this section we will apply directly the results of the 
previous sections to study establishment of entanglement 
over large scales in lattices with hierarchical geometry. 
These are lattices that iterate certain geometric struc- 
tures, so that at each level of iteration the number of 
nodes or the number of neighbors changes. Unfortu- 
nately, we do not know how to find optimal strategies 
for such lattices; we restrict our considerations to show 
that one can establish perfect entanglement in a finite 
number of steps at some iteration level. This perfect en- 
tanglement can be swapped further to the lowest levels 
of iteration, i.e. to the largest scales, which can be con- 
sidered as the largest geometrical distances. One should 
stress that classical percolation strategies work for hier- 
archical lattices as well as for 2D or 3D ones. Neverthe- 
less, every percolation strategy relies on conversion of all 
bonds to singlets with a certain probability p and then es- 
tablishing a perfect entanglement between two nodes on 
the large scale with a probability 9^{p). The latter for- 
mula expresses the fact that both nodes have to belong to 
the percolating cluster, which happens for each of them 
independently with probability 9{p) This probabil- 
ity is always smaller than one, except for the trivial case 
of p = 1. 



A. "Diamond" Lattice 

We start considering the so-called "diamond" lattice, 
which is obtained by iterating the operation presented in 
Fig. [71 in which a single bond (two qubits and one entan- 
gled state) is replaced by four bonds forming a diamond 
shape (four pairs of qubits and four entangled states). 
We prove that for sufficiently large initial entanglement, 
one can establish perfect entanglement on large scales 
{i.e. on some lower levels of iteration) in finite number of 
steps. 

We assume that the lattice is formed by very many 
iterations, and that all bonds correspond to entangled 
states \if) — ipo |00) -I- (pi |11). Our aim is to perform 
measurements in a recursive way and demonstrate the 
for sufficiently small ipo it is possible to establish per- 
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FIG. 7: The diamond lattice is formed by iterating the follow- 
ing operation: a single bond (two qubits and one entangled 
state) is replaced by four bonds forming a diamond shape 
(four pairs of qubits and four entangled states). After K iter- 
ations, the nodes A, B, C, D have 2^ links, the nodes on the 
next level 2^^^ links, etc. 



feet entanglement on the lowest level of the iteration hi- 
erarchy, i.e. between the "parent" nodes A and B. In 
order to keep the form of the network unchanged dur- 
ing the recursive measurement we will apply the WCE 
strategy to the nodes analogue to C and D, staring from 
the highest (last) iteration level. After applying WCE 
we obtain with probability 1 a pair of entangled states 

I^A) = V'o |00) + Vi 111), with ^0 = (1 + \/l - 16^o'^D/2- 
This pair can then be distilled with probability 1 to a 
new two-qubit entangled state \^'), see Eq. (P^ : 

^'q = max |i, i (l + y^l - l&^l^t) I ■ 

Denoting now the SCP by = 2(pi, we rewrite the re- 
cursion as 

E' = 2(1 - V'o) = 1 + (2 - - Vl - (2 - EfE"^. 

This recursion (see Fig. [8]) has one nontrivial unstable 
fixed point Eth, and two trivial stable fixed points E ~ Q 
and E ~ \. The latter is achieved in a finite number of 
steps provided the initial E > Eth ~ 0.349. Note that 
Etii is strictly smaller than E* = 2(1 — (^3) from Eq. ((30)) . 
For E > E* , E' is equal to 1, i.e. the singlet is achieved 
in one step. 

B. "Tree" Lattice 

Similar results hold for the simplest possible "tree" 
lattice: a double Cayley tree lattice with branching factor 
two (see Fig. [3^). 

Let us denote the initial SCP of all bonds by Eq . The 
strategy is depicted in Fig. [HI the nodes in the middle of 
the tree perform the WCE. This prepares two two-qubit 
states between the neighboring nodes with the entangle- 
ment, measured by the SCP, equal to Ei. These two 
states are then converted with probability 1 into a two- 
qubit state with E = min{l, 2(1 - (1 - Ei/2)^)}, which 
will undergo recursive transformations (Fig.Oj). We per- 
form then the WCE on one of the 3 connected bonds, and 
obtain Ei ^ 1 - ^1 - Eo{2 - Eq)E{2 - E). Then, the 
WCE is applied to the remaining pair of bonds yield- 
ing £'n = 1 - ^1 - Ea{2 - Eq)Ei{2 - Ei). Finally, the 




0.2 0.4 0.6 0.8 1 



Initial entanglement E 

FIG. 8: Recursion relating E on the higher level of lattice 
iteration to E' at the lower level of iteration in the diamond 
lattice. Each iteration consists of the following steps: (i) WCE 
and (ii) the two resulting two-qubit states are transformed 
with probability one into a two-qubit state of the same SCP. 




FIG. 9: (a) Tree configuration; (b) The nodes in the middle 
perform WCE. This creates two two-qubit states between the 
neighboring nodes. These states are transformed into a two- 
qubit state of the same SCP. The process is iterated until a 
perfect singlet is established between the two ends of the tree. 



optimal singlet conversion is applied to the pair of En 
bonds obtained from the two different but neighbouring 
branches of the tree, yielding E' = min{l, 2(1 — (1 — 
i?ii/2)^)}. The recursion relations can be rewritten as 

E' = F{E;Eq). (33) 

This recursion depends explicitely on Eq. It is easy to 
see that since the WCE does not increase the SCP, the 
recursion (j33p will have: i) only one trivial stable fixed 
point E = if Eq < Eth, ii) three fixed points otherwise: 
stable 0, unstable E and stable 1 otherwise. Here if we 
start with E > Eth we will end up in one step with 
E' = 1. The threshold value is obtained by solvi ng the 
equation 1 = 2(1 - (1 - (1 - ^^1 - E^{2 - Eoy)/2y), 
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and is given by 



(a) 



«0 — CL 
O O 



/^O — CL 
O O 



th 



1- \/l- V2(\/2-l) wO.7. 



(34) 



6^4 4^6 



VI. GENUINE 2D LATTICES 

In this section we consider genuine 2D lattices when 
the number of nodes is big. On the one hand, we apply 
the methods and observations of the previous sections to 
these large lattices. On the other hand, we reconsider 
the various variants of the methods employing classical 
and quantum percolation techniques. 



A. "Centipede" in Square Lattice 

As another example of the power of recursive mea- 
surement methods of the previous section, we consider 
a wide strip of a 2D square lattice and the "centipede" 
figure within it (see Fig. [TUb). Let us denote the ini- 
tial entanglement as E^), and the entanglement at the 
end bond of a "leg" by E. We then apply the following 
measurement scheme to the ends of each of the legs of 
the centipede, see also Fig. [TUl i) We apply the WCE to 
Eq and i?, replacing these two bonds by one with Ei = 
1 - a/1 - Eo(2 - Eo)E(2 - E); (ii) we repeat it with the 
other vertical bond obtaining thus a pair of states at the 
horizontal end of the leg: one with entanglement Eq and 
the other with Eu = 1 - ^\ - Eq{2 - Eq)Ei{2 - Ei); 
(iii) the resulting pair is then distilled with probability I 
to a two-qubit state with £" = min{l, 2(1— (1— £'o/2)(l — 
i?ii/2))} — F{E; Eq). This situation is somewhat similar 
to the case of the tree lattice from the previous section, 
but not completely. The recurrence relation depends ex- 
plicitly on Eq and has always a nontrivial stable fixed 
point E > Eq. This fixed point, however, is strictly 
smaller than 1, when Eq is small. In the first case, al- 
though we do concentrate more entanglement along the 
"spine" we still have to face the problem that the spine 
is a ID network, and will exhibit exponential decrease of 
probability of establishing the perfect entanglement for 
large distances fTTl. On the other hand, the stable fixed 
point is simply E = I, provided Eq is large enough. In 
this case a perfect singlet is achieved in a finite number of 
steps, and the singlets from all legs can be concentrated 
at the spine of the centipede with probability 1. Obvi- 
ously, all that implies that the width of the strip of the 
2D lattice (equal to twice the length of the centipede leg) 
can be finite: it must be just larger than the number of 
steps necessary to get a perfect singlet. 

The condition for the threshold value of Eq is easy to 
derive: we have to put E = lin the above recurrence and 
solve 1 = 2(1 - (1 - Eth/2){1 + ^/T 
which gives Eth ~ 0.649. 



EU2 



6 6 
.0 — o" 



A A 
P — o. 



— CP3 — o"o — o"o — o"o — o"o — o"o — oxo — o 



TP — cr 
o o 



% — cr 
o o 



! I ! 



% — cP 
o o 



w w w w 



(b) 



f 


E 






Eo 






En 








5f 



Eo 



Eo 




FIG. 10: (a) "Centipede" with its "legs" and "spine"; (b) 
recursive measurement scheme; note that the method can be 
equally well applied also in higher dimensions. 
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FIG. 11: Recursion relation for the centipede lattice. Only 
when the entanglement Eo is larger than a threshold Eth, a 
trivial stable point at iJ = 1 appears. 



B. Percolation Strategies 

In Ref. [m we have pointed out that one possible strat- 
egy for entanglement distribution is to convert locally 
all bonds with probability p into singlets and then per- 
form entanglement swapping. This strategy can then be 
linked to classical percolation theory, so all the known 
results of this field can be applied to our quantum sce- 
nario. The natural question is whether the thresholds 
defined by classical percolation theory are optimal or en- 
tanglement percolation represents a related but different 
theoretical problem where new bounds have to be ob- 
tained. This of course equivalent to determine whether 
the measurement strategy based on SCP is optimal in the 
asymptotic regime. Here, we construct several examples 
that go beyond the classical percolation picture, prov- 
ing that the classical entanglement percolation strategy 
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lattice 


Pc 


triangular 




= 2 sin(7r/18) fa 0.347 


square 




= 0.5 


honeycomb 


P? 


= 1-2 sin(7r/18) ^ 0.653 



TABLE IV: Classical (bond) percolation thresholds pc for 
some regular lattices. 



is not optimal. The key ingredient for the construction 
of these examples is the measurement strategy previously 
obtained for the ID one-repeater configuration that max- 
imizes the SCP. 



1. Honeycomb Lattice with Doubled Bonds 

The first example (already discussed in Ref. [Tlj ') 
considers a honeycomb lattice where each node is 
connected by two copies of the same two-qubit state 
|(^), see Fig. [Hi. 

The simplest strategy consists in trying to convert all 
bonds of the doubled honeycomb lattice into singlets, 
and then applying the classical entanglement percola- 
tion. The percolation threshold of this doubled lattice 
is not difficult to calculate: at the critical point, the 
probability that, at each edge, at least one conversion 
is successful has to be equal to the percolation thresh- 
old of the simple honeycomb lattice; if both conversions 
are successful we simply discard one pair. We thus have 
pP = 1 — (1 — 2(pi)^, hence the percolation threshold is 



(see Tab. |TV 



2(^1 = 1 — W 2 sin 



0.411. 



(35) 



We define the classical entanglement percolation strat- 
egy as i) converting in the best possible way all bonds 
shared by two parties into one singlet and ii) applying the 
entanglement percolation. If, as above, the Schmidt co- 
efficients of the two-qubit state are fo> (pi, the SCP of 
is given by pok = 2(1 — (/Jg). We choose this conver- 
sion probability to be equal to the percolation threshold 
for the honeycomb lattice and get 




— « 0.358. 



(36) 



We now show that another strategy yields a better per- 
colation threshold: some of the nodes, see Fig. [T2k . per- 
form the optimal strategy for the SCP, mapping the 
honeycomb lattice into a triangular lattice, as shown in 
Fig. [T2b . What is important is that the SCP for the new 
bonds is exactly the same as for the initial state |(^), that 
is 2(^1. We choose it to be equal to p^, so that 




2 sin 



0.347, 



(37) 



FIG. 12: Each node is connected by two copies of the same 
two-qubit state The nodes marked in (a) perform the 

measurement optimal according to the SCP. A triangular lat- 
tice (b) is obtained where the SCP is the same as for the state 
|(^). Classical entanglement percolation is now possible in the 
new lattice. 



which proves that the classical entanglement strategy is 
not the best one. 



2. Asymmetric Triangular Lattice 

The second type of examples, although less symmetric, 
is generic and has a totally different character than the 
previous one. For simplicity, we show the argument in the 
case of a triangular lattice, but the same reasoning can be 
applied to other geometries. Consider the triangular lat- 
tice of Fig. 113b . Solid lines correspond to two-qubit pure 
states \tf) while dashed lines correspond to states |<^) 
that are less entangled, i.e. cpo > ipo. We choose the first 
state such that Pok = '2(pi satisfies p^ < Pok < \/p^- If 
= \ip), the classical entanglement percolation strat- 
egy works. However, we choose this second less entan- 
gled state such that its SCP is small enough to make 
the classical entanglement percolation impossible. This 
state always exists. Indeed, note that when ipi — > 0, these 
states can simply be removed from the lattice, and clas- 
sical entanglement percolation fails because of p^j. < p^ . 
It is now rather straightforward to construct a success- 
ful entanglement percolation strategy: the state \(p) are 
discarded and the optimal strategy for the one-repeater 
configuration and the SCP is performed. The lattice is 
then mapped into a new triangular lattice keeping the 
conversion probability of the first, more entangled, state, 
see Fig. 113b. Classical entanglement percolation can now 
be applied to this new lattice, since pok > Pc- 



C. Doubling the Square Lattice 

The final example deals with a square lattice and has 
yet another character. Here we replace every second pair 
of horizontal bond by a single one using the optimal SPC 
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(a) (b) 




FIG. 13: The triangular lattice consists of two different en- 
tangled states \ip) and |(^) for the solid and dashed lines, re- 
spectively. The less entangled states, \(p), are discarded and 
some of the nodes perform the optimal measurement accord- 
ing to the SCP. A new triangular lattice is obtained, governed 
by the SCP of 




FIG. 14: (a) Measurements necessary to double the square 
lattice: the marked nodes apply the optimal one-repeater 
transformation along the vertical and horizontal directions, 
(b) Resulting pairs of disjoint square lattices with lattice con- 
stant doubled; we want to establish perfect entanglement be- 
tween any two neighboring points A, A' versus B and B' . 
A and A' [B and B') are neighbors but belong to different 
lattices. 



strategy, which as we know from Sec. |TTT] does not change 
the SCP on average, replacing, however, pure states by a 
known mixture. The same is done with every second pair 
of vertical bonds. In effect we replace the original square 
lattice by two disjoint lattices with the lattice constant 
twice bigger than the original one, but the same SCP (see 
Fig. [H]). Now we consider the following problem: we are 
interested in establishing entanglement between any of 
the two neighboring nodes A, A' and B, and B' at large 
distances. 

In the case of doubled lattices the calculation is sim- 
ple: the pairs {A,B) and {A' ,B') belong to two disjoint 
lattices, and the probability that, say, (A, B) belongs to 
the percolating cluster is equal asymptotically to d^{p). 
The probability that at least one of the pairs belongs to 
the percolating cluster is thus 

PdouUc = 2 9^-9^ = e^ (2-9^). (38) 

This probability has to be compared with the probability 
that at least one of the pairs (A, _B), {A',B), {A,B') or 



{A' , B') belongs to the percolating cluster in the original 
square lattice. The latter probability is asymptotically tt^ 
where tt is the probability that A or A' (or equivalently 
B or B') belongs to the percolating cluster C. Thus we 
have 

TT = P[Aot:A' e C] 
= P[A e C] + P[A' e C] - P[A, A' e C]. 

In order to estimate the last term in the above expression, 
we use the FKG inequality ([1^, see also i_24^|). To state 
it, we first define an event described in terms of a percola- 
tion configuration to be increasing if it has the property 
that, once it holds for a certain bond configuration, it 
holds for all configurations obtained by adding bonds to 
the initial one. FKG inequality says that any two such 
events are positively correlated. The events {A e C} and 
{Ao-oA'} ("A and A' are connected by a path of maxi- 
mally entangled bonds" ) are clearly increasing and, since 
their intersection is the event {A, A' e C}, it follows that 

P[A,A' eC]>P[AeC]P[Ac^A']. (39) 

Denoting P[A^A'] by r, we thus have tt^ < 9'^ 
Therefore, doubling the square lattice is a better strategy 
than the classical percolation, i.e. tt^ < Pdoubic, whenever 

(2-r)2 < 2-6*2. (-4Q<) 

We believe that this inequality is fulfilled for all p, and 
we show it when p is just above the percolation threshold 
p° = 0.5, so that 9 tends to zero. To this aim we have to 
demonstrate that 2 — r < ^/2. We may try to estimate 
r from below by considering the six shortest trajectories 
connecting A and A': the most direct ones, two 2-edge 
paths, and the two pairs of 4-edge paths around the ad- 
jacent squares. One finds 

T > 2 (p2 + 2/ - 2p5) - (p2 + 2/ - 2p5)2. 

Unfortunately, for p — p^ this estimate is too small, since 
it gives only 2 — r < 1.473 . . . and adding further paths 
becomes then technically tedious. We therefore turn to 
the standard numerical Monte Carlo method, generating 
the shortest paths automatically while the longer ones are 
generated using the Monte Carlo sampling. For p > p° 
the convergence is exponential: if we plot a subsequent 
estimate of r as a function of the maximum cluster size 
allowed in the Monte Carlo sampling, it approaches the 
final value exponentially fast for large clusters. As ex- 
pected, the convergence is algebraic at p = p°: the esti- 
mate of T approaches its final value as a power of the clus- 
ter size. A power law fit and a comparison with the val- 
ues just above the percolation threshold give with a very 
good accuracy r ~ 0.687 and hence 2 - t ~ 1.313 < \/2, 
q.e.d. 

This is yet another result which does not have a clas- 
sical analogue, showing how quantum mechanical mea- 
surements allow to increase percolation probability. 
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VII. CONCLUSION 

In this paper, we have considered the problem of en- 
tanglement percolation through pure-state quantum net- 
works. We have first focused our investigations on small 
quantum networks. Even for these particularly simple 
systems, interesting and unexpected properties have been 
pointed out. One of the main result is the description of a 
Bell measurement by its outcome probabilities only (Re- 
sult [T|) . This has allowed us to maximize the different 
figures of merit introduced at the beginning of the pa- 
per. We have shown, then, that Bell measurements do 
not yield in general the optimal protocol, even for a chain 
consisting of only two repeaters. 

The results for small lattices have later be used as 
building blocks for entanglement percolation protocols in 
asymptotically large lattices. We have provided several 
examples illustrating some of the properties character- 
izing these lattices: recursive relations, classical entan- 
glement percolation examples of lattices were quantum 
effects allow going beyond classical percolation. 

In general, little is still known about the problem of 
entanglement percolation, that is, the distribution of en- 
tanglement through quantum networks. In the pure-state 
case, it would be interesting to derive lower bounds to the 
amount of entanglement between the nodes such that en- 
tanglement percolation is possible. The main question, 
however, is to extend these results to the mixed-state 
scenario, providing examples of entanglement percolation 
protocols for lattice with mixed-state bonds. 
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that Pm = Pmin + Pmax (1 " km) with km & [0, 1] . One 

of the four equations of the system {pm} = {xm} will be 
dependent of the other three: if we can find three orthog- 
onal vectors such that pm = Xm for, say, m = 1, 2, 3, 
then the fourth one is fixed (up to a sign) with, obviously, 
Pi — X4^. Let us write these three states as 

Mm = {^/k^ cos(0„i)' V^m sin(6'm), 

yjl-km cos(tj„i), \/l-km sin(tj„i)), 

where km = (Pmax-2^m)/(Pmax-Pmin)- By coustructiou, 
these vectors are normalized and satisfy pm = Xm- We 
now have to prove that there always exist some angles 
Om and LUm such that these three vectors are orthogonal. 
Without loss of generality we order the /c's such that 
^ > ki > k2 > k^ > ki> 0. Since the probabilities add 
up to 1 and that Pmin +Pmax = 0.5 we have 

/ci + /S2 + ^3 + ^4 = 2. (A.l) 

Introducing the notations k',^ = 1 — km^ da = 01 — 02, 
9b = 01 — 03, uJa = uJi — UJ2, uJb = 1^1 — W3 and using 
the identity cos(a:) cos(?/) -I- sin(a;) sin(y) — cos{x — y) , the 
conditions of orthogonality read 

{0 = cos(6'q) -f ^Jk'^k'^ cos(wq) 

= y/kiks cos{0b) + \/k[k'^ cos{LUb) 
= Vhh cos{0a - 0b) + y^jps cos(a;a - cut) 

(A.2) 

The cases km = or km — 1 for some m can be triv- 
ially solved, so we consider fc^ 7^ and k'^ ^ 0. We 
have four parameters 0a,b and oja.b which can be freely 
chosen in the interval [0, tt], but the two inequalities 
Vkik2 > y/k[k^ and ^/kik^ > \/k[k^ impose the con- 
straints 0a e [0'a, TT - 6**] and 0b e [01, TT - 01], with 
0*a b e [0, f ] such that ^/k[k^ cos(6'*) = ^/cffc^ and 
"v/fcifcs cos{0l) = y^k[k^. Thus cos(wa— Wf)) & [—1, 1] and 
cos{0a - 0b) e [~cos(6'* -I- 01), 1]. Then, one can verify 
that there always exists at least one solution of Eq. (jA.2[) . 
except when —^^2^3 cos(0* + 0^) > \/kl^^, what never 
happens. In fact, suppose that this last inequality holds 
and rewrite it in terms of k\ , k2 and k^ only. After some 
tedious algebra and using some trigonometric identities, 
one finds that the inequality fci + fc2 + fca > 2 holds, but 
this is in contradiction with Eq. (jA.ip . which concludes 
the proof. □ 



APPENDIX A: PROOF OF RESULT [T] 

We prove here that there always exists a Bell mea- 
surement which yields outcome probabilities Pm equal 
to Xm, when these values add up to one and lie in the 

interval [pmin, Pmax]- 

Proof. (By contradiction). Let us write {fim} the four 
states of the Bell measurement in the magic basis. Be- 
cause the matrix {p-m.i) is orthogonal, the conditions on 
X are clearly necessary. In fact, we know from Eq. ([S]) 



APPENDIX B: SCP OF ZZ MEASUREMENTS 
ON A ID CHAIN 

Even if the number of outcomes grows exponentially 
with the number of repeaters, one can keep track of all of 
them in an efficient way. In fact, after any number n < N 
of entanglement swappings in the ZZ basis, any possible 
resulting state has the form (up to local unitaries) 

|m) = ^=L=(y^|00) + v/;^|ll)), men. 
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n = n = 1 



FIG. 15: "Tree view" of the labels m and their corresponding 
probabilities after the first two measurements. For symmetry, 
we choose the root of this tree corresponding to n = — 1. 



We prove this by induction on n, the case ti = corre- 
sponding to the initial state m = 1. Suppose that the 
result holds and that we got the state |m) after n < N 
measurements. It is easy to show ft-om Eq. (fT8|) that an 
entanglement swapping in the ZZ basis on |m) (g) \ip) is 
described by 

I |m + l) with probability p+ (m) /t^,n 
\m) 1-^ < (Jj-l) 
I I |to — 1| ) with probability p_(m), 

where p+{m) = + ip'^+^)/{ipl^ + ipl') and 

p^(m) ~ 1 — p_|_(rn), q.e.d. 

The first step to calculate the SCP of this protocol is 
to compute its variation after a ZZ measurement. Con- 
sidering that the set {(pi, \mi)), i = 1, . . . describes 



all the resulting states of n measurements, and writing 
A±(m) the smallest Schmidt coefficient of |m ± 1), the 
new SCP reads 

I 

1=1 

= 4'^)-((/.o-^i)p(m = 0,n), (B.2) 

where p(m = 0, n) stands for the probability of getting 
the state |m = 0) after n measurements. Since this prob- 
ability is not zero for n odd only, it results that the SCP 
decreases for n even only. We have now to calculate the 
probability p{m = 0, n) of getting a singlet after n mea- 
surements: it is the weighted sum over all possible paths 
r that go from the root node m = to the node m = 
at position n in the tree drawn in Fig. 1151 We notice 
that the weight, w, depends on n only and not on F. 
This is indeed the fact since p^{m) p-{m + 1) = Lp^ ip\ 
for all m and because we have to go up in the tree as 
many times as we have to go down. Thus, for n odd we 
have — (ipQi^xf-'^^^'^^'^ and using basic combinato- 

rial analysis one finds that p(m = 0,n) = (ipQipx)^ (^^^ , 
with fc = i(n + 1) e N. Finally, denoting by \x\ the in- 
teger part of x, the general expression of the SCP for a 
chain of N repeaters reads 

4^' = 1 - (¥^0 - ^i) E ('^o^i)' J- (B-3) 
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